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Abstract. We investigate the coupling between a quantized electromagnetic field in 
a cavity resonator and a Coulomb interacting electronic system in a nanostructure in 
an external magnetic field. Effects caused by the geometry of the electronic system 
and the polarization of the electromagnetic field are explicitly taken into account. 
Our numerical results demonstrate that the two-level system approximation and the 
Jaynes-Cummings model remain valid in the weak electron-photon coupling regime, 
while the quadratic vector potential in the diamagnetic part of the charge current 
leads to significant correction to the energy spectrum in the strong coupling regime. 
Furthermore, we find that a coupling to a strong cavity photon mode polarizes the 
charge distribution of the system requiring a large basis of single-electron eigenstates 
to be included in the model. 
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1. Introduction 

In the last decade there has been increasing interest in systems capable of generating 
quantized fields containing a preset number of photons. Manipulation of the state 
of scalable light-matter coupled quantum systems is one of the key issues for their 
implementation for optomechanical systems [U |2] or quantum information processing 
devices [31 H] . However, searching clear evidence of light-matter coupling nonlinearity 
is still a challenge. To this end, one has to reach a strong light-matter coupling regime 
for optically driven systems in high quality micro- cavities [51 |6], and demonstrate its 
single-photon characteristics [TJ Ej. Flexible experimental design of circuit quantum 
electrodynamics offers great potential for practical device applications to explore strong 
light-matter coupling at microwave frequencies [91 [101 [TTJ [121 02] • 

To describe the interaction between matter and the photons of a quantized 
electromagnetic field, the Jaynes-Cummings (JC) model is often applied [H]. The 
JC model describes the interaction between a two-level system (TLS) and a single 
field mode. It is a fundamental model in quantum optics and quantum information 
science [15J. For a TLS with energy level spacing AE, coupled with strength £jc to 
a resonator with photon energy hut, the JC model is valid when both the detuning 
5 = \ftw — AE\ is sufficiently small and the light-matter coupling strength is much 
smaller than the photon energy (£jc <C hoj). The dynamics can then be obtained by 
the JC model [16] and the energy spectrum can be solved exactly if the rotating wave 
approximation (RWA) is applied [TT] . 

Single modes of the electromagnetic field can be treated as the population of a 
field oscillator with different Fock states (states with certain number of photons). It 
was discovered that a three-level problem, called the coupled-channel cavity quantum 
electrodynamics model [18] . can be exactly transformed to a two- level one for arbitrary 
detuning [19] , in which the eigenstates of energy and orbital angular momentum can be 
explicitly expressed in terms of the Fock states. 

More recently, utilization of the giant dipole moments of intersubband transitions 
in quantum wells [201 121] has enabled researchers to reach the ultrastrong light-matter 
coupling regime [221 1231 121] ■ In this regime, the JC model is not applicable and the 
coupling mechanism has to be explored beyond the JC model [251 EEl E7] . Despite the 
above mentioned experiments, a study of the coupling between electrons and cavity 
photons with a specified nanostructure geometry in a perpendicular magnetic field is 
still lacking. 

In this work, we investigate the interplay of the dynamics of correlated electrons 
in a nanostructure to the quantum field of a rectangular cavity resonator subject to an 
external magnetic field. By performing numerical computations we demonstrate how 
the electron-photon coupling influences an electronic system embedded in a quantized 
photon field. The TLS approximation and the JC model will be examined in both the 
weak and the strong coupling regimes as well as the effects of the diamagnetic part of 
the charge current in the electron-photon interaction term, which the JC model lacks. 
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2. Model and Theory 

The system under investigation is a two-dimensional electronic nanostructure exposed 
to a quantized electromagnetic field of a cavity resonator and a static (classical) external 
magnetic field at a low temperature. The electron-photon coupled system can be 
described by the many-body Hamiltonian 

H = H c + H EM + H C _ EM , (1) 

where the first term describes the electronic system including the magnetic-field modified 
kinetic term H® and a Coulomb interaction term i?coui> namely 

H e = H® + Hcoui- (2) 

The second term Hem in Eq. (OQ) represents the electromagnetic field in a cavity resonator 
and the third term H e _ EM contains the coupling between the electronic system and the 
quantized electromagnetic field. 

The electronic nanostructure is assumed to be fabricated by split-gate configuration 
in the y-direction, forming a parabolic confinement with the characteristic frequency Qq 
on top of a semiconductor heterostructure. The ends of the nanostructure in the x- 
direction at x — ±L x /2 are etched, forming a hard- wall confinement of length L x . 
Thereby, a closed electronic narrow constriction is created in the 2D electron gas. 
The external classical magnetic field is given by B = Bz with a vector potential 
A cxt = (—By,0). Hence, H® can be expressed in second quantization as 

H *= /^ t { 2 ^ + ^ m ^ 2 }^' (3) 

where 7r = p + eA ext /c is the mechanical momentum, m* is the effective electron mass, 
and 

V = X;V<(r)d i , ^ = J>*( r )4 ( 4 ) 

i i 

are fermionic field operators, with di the annihilation- and d\ the creation operator for 
an electron in the single-electron state \i) corresponding to ipi. Through out this paper 
we use Latin indices to label single-electron states (SESs) as well as Hilbert state vectors 
and Greek subscripts to label many-electron Fock states (MESs). Rationalized Gaussian 
units are used exclusively and e denotes the positive elementary charge. In Eq. (j3J), ^i(r) 
can be any complete orthonormal set of functions with the correct boundary conditions. 
However if we use the SESs of H® as a basis, H® is diagonal and simplifies to 

H° e =Yl E i d tdi, (5) 

i 

where is the energy of the SES i, associated to the eigenfunction ^(r). The SESs 
are computed numerically in a functional basis using a straightforward diagonalization 
method. 



Quantum magneto- electrodynamics of electrons embedded in a photon cavity 



4 



We can write the Coulomb interaction term in the second quantized form 
1 
2 



ffcoui = -z ^2 I ^coui I rs) d\d)d s d r , (6) 



Vco U i(r,r)= , , (7) 



ijrs 

where the Coulomb interaction potential can be written as 

e 2 /K 
|r — r'| + rf 

with k denoting the relative dielectric constant of the material and r\ an infinitesimal 
convergence parameter. The Coulomb matrix elements in Eq. ([6]) are thus expanded in 
the basis of the SESs involving the integration with respect to the observing location r 

(u'|Vboui|rs) = J dr V;(r)X ir (r)V,(r) (8) 

and the integration with respect to the source location r' 

X jr (r) = J drV^rO^Coui^OVv^') . (9) 

We use the SESs to construct many electron Fock states which obey H®\n) = 
Ey\y) and employ an exact numerical diagonalization method to obtain, in the 
(noninteracting) basis the Coulomb interacting eigenvectors \fi) which satisfy 

(H® + .ffcoui)!/-*) = E^\fi). We write the Coulomb interacting eigenvectors in terms of a 
unitary transformation [291 SO] 

Im) = X)H»> (io) 

u 

which is obtained in the diagonalization process. Then, in the new interacting MES 
basis {|/x)}, we can rewrite the full electron Hamiltonian as 

H e = J2\^)EM- (11) 

The cavity electromagnetic field is described by the quantized vector potential A in the 
radiation (Coulomb) gauge. The free field Hamiltonian is simply 

H EM = huxra , (12) 

where u is the frequency of the resonant cavity mode, a) and a being the creation and 
annihilation operators for photons. The last term in Eq. (JTJ) describing the interaction 
between the electrons and the quantized electromagnetic field is given by 

1 

c J '" 2m*c 2 
in which the charge current density is defined by 

j e = {^W) + (ir'VW ' ( 14 ) 

and the charge density p e = —eip^ip. In this paper we will show that the A 2 term in 
Eq. ( !T3|) significantly affects dynamical features in the strong electron-photon coupling 
regime. The electronic nanostructure is placed in a rectangular cavity forming an 



#g-em = ~- I j c ■ A - — - - / dv p e A 2 , (13) 
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electromagnetic oscillator with hard-wall boundaries at —a c /2 < x, y < a c /2 and 
—d c /2 < z < d c /2 with cavity volume V c = a 2 c d c . The proposed electromagnetic 
oscillator is a single planar rectangular cavity with mutually locked dual anti-phase 
outputs, in which the electronic nanostructure is oriented in the z = plane with the 
center at (x,y) = (0,0). In the following, we will consider only transverse electric 
(TE) modes (E z = 0), where the electric field E is perpendicular to the direction of 
propagation. The cavity supplies a monochromatic wave stabilized in the TE n mode 
with longitudinally polarized electric field along x-direction, or in the TEioi mode with 
transversely polarized electric field along y-direction. In the Coulomb gauge, the vector 
potential of the electromagnetic field takes the form 



A(r) 






' 7TX \ 




cos 








r< 




COS 








K a c J 





COS 



^)(o + ot) (15) 



with the upper component denoting the TE n mode and the lower one representing the 
TEioi mode. The Cartesian unit vectors are e x and e y . 

We assume that the size of the cavity is much larger than that of the nanostructure, 
that is L x <C a c , d c . Utilizing this condition we can approximate the cosines in Eq. ( fT5l) 
by unity and take A outside the integrals in Eq. (fT3|) and obtain 

i? e _EM = £ c ^2 d i d i 9ij ( a + flt ) ( 16 ) 

l\ 1 , : 

- ya'a' + aa) 



si 

MX, 



L w 



a) a 



2 7 2 



where Vt w = a/^o + w c i s ^ ne effective confinement frequency, u c = eB/(m*c) the 2D 
cyclotron frequency, S c = eAa w Q w /c the characteristic energy for the coupling between 
electrons and cavity photons, a w = \/h/ (m*Q w ) the characteristic oscillator length scale, 
and M e = Yli d\di the electron number operator. In Eq. f fT6|) we use the dimensionless 
coupling gij between the electrons and the cavity modes defined by 

9H = ^J dr t^W {(e • Tr)^-(r)} 

+ {(e- 7 r)^(r)}> J (r)], (17) 



with e ■ 7r = e x ii x + e y 7i y . The first and the second terms in Eq. (1161) contribute, 
respectively, to the linear and nonlinear optical excitation energy spectra, which will be 
explored later. 

In our theoretical consideration, we will formally consider all the higher-lying 
photonic modes, truncating the infinite matrix in order to retain enough modes to 
reach sufficient convergence. To obtain convergent energy spectrum from the total 
Hamiltonian Eq. (pQ), all the resonant and antiresonant terms in the photon creation 
and annihilation operators will be taken into account with arbitrary detuning. In 
the case of resonance, the condition of the vacuum Rabi frequency to the cyclotron 
frequency ratio larger than one implies that the higher photonic modes are coupled to 
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the transition, and the diamagnetic A 2 term of the electron-photon coupling in Eq. (115]) 
becomes dominant [28J. 

In the MES basis {|/i)}, we can rewrite the electron-photon interaction Hamiltonian 

as 



H e -EM = £ c \y){y\V^d\djV\v){v\ g {i {a + a)} 



El . 

v\ 



w pvj 



x 



o)a + -J + - [aa + aW) 



The energy spectrum of the total Hamiltonian involving both the electron-photon 
and electron-electron interactions has to be obtained from the many-body space of 
the interacting electrons {|/i)} and the Fock-space of photons {|M p h)}, namely the 
uncoupled electron-photon many-body states (MBSs) \p) = <S> \M p t). Using 
the noninteracting electron-photon MBSs, performing diagonalization of the total 
Hamiltonian ([TJ), the interacting electron-photon MBSs can be expressed as 

V 

It is important to note that, in arriving to Eq. (fl9l) we have performed truncations 
to a basis two times. The first one is when we only use the iVsES single electron 
states to construct our Fock-space basis The second is when we only use N MES 

Coulomb interacting states and Nem photon states to construct the joint photon- 
electron many-body basis <g> |M ph )}. A third truncation is likely to be needed if 
one needs to apply operations involving W many times, for example in time dependent 
calculations [3]]. However in this paper, only static properties are calculated such as 
the energy spectrum and charge densities of the many-body eigenstates. Therefore a 
third truncation is not necessary. 

The electron charge density operator Q(r) in the electron-photon coupled system 
in the second quantized form is 

Q( r ) = -eJ2^(r)MMd r (20) 

By taking trace of the operators in the coupled MBS space {|/t)}, we obtain the 
correlated many-body charge distribution (Q(t)) for the electrons in the nanostructure 
interacting with the photon field 

(Q(T))=Tr{p(t)W*Q(r)W}. (21) 

It should be mentioned that the total density operator p(t) contains information of the 
interacting many-electron system and the monochromatic photon modes as well. 

Below we shall demonstrate our numerical results displaying tunable dynamical 
interplay features between the interacting electrons and the quantized photon field with 
either x or the y polarized electric components in an external magnetic field. We mention 
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in passing that in the construction of the JC model we assume only two electron states 
to be active in the system. However, when we take a realistic geometry into account, a 
sufficient number of electron states will be used for guaranteeing numerical convergence. 

3. Results and Discussion 

In order to explore the dynamical features of an electronic system coupled to a single- 
mode quantum photon field, a simple TLS model is often employed [32j [331 EH] ■ The 
energy spectrum of the TLS can be used to employ the JC model [TS] . Although the 
JC model has been applied beyond the RWA, a rigorous analysis for the validity of 
the JC model in a realistic system remains unexplored. It is thus useful to consider 
an electronic nanostructure coupled to a cavity photons taking into account realistic 
geometrical effects, by performing numerical computation beyond the JC model for 
comparison. 

We assume that the electric nanostructure is fabricated by GaAs-based 
semiconducting materials with electron effective mass m* = 0.067m e and the background 
relative dielectric constant k = 12.4. The electronic nanostructure is modeled as an 
infinite square potential well of length L x = 300 nm along the x-axis and transversely 
constricted by split-gates with a parabolic confining strength h~Q = 1.0 meV along the 
y-axis. A uniform and static magnetic field B = 0.1 T is applied along the z-axis. In 
order to obtain the energy spectra of the closed electron-photon system, we have used 
^Vmes — 200 electron states and N EM = 20 photon states. The two electron MESs have 
been computed with iVsEs = 50 and the single electron states with N$es = ^mes- No 
calculations for three or more electrons are needed in this paper. 

The MBS energy spectra of the electron-photon states are shown in Figs. [T] and 
[2] for the case of x and y-polarizations respectively, corresponding to the TE n and 
the TEioi modes. The single photon energy is hw = 0.4 meV, hence the M p h photon 
states with no electrons have the energy = Mp^hu. The horizontal axis £ c denotes 
the strength of coupling between the electrons and the photons. In the absence of the 
electron-photon coupling S c = 0, both the cases of x- and ^/-polarizations manifest the 
same energy spectra. 

In the weak electron-photon coupling regime S c < O.lhco, the energy spectra of a 
linear approximation neglecting the A 2 term (red) are only slightly shifted from the full 
numerical results including the A 2 term (blue), and the pure photon states (green) retain 
the same energy. However, the linear approximation neglecting the A 2 term becomes 
inaccurate when the coupling strength S c is comparable to the driving photon energy 
hu and the characteristic energy level spacing of the electronic system. 

Our full numerical results clearly show that the lowest two many-body energy states 
at around E = 0.56 and 0.75 meV are not sensitive to the polarization of the quantized 
electromagnetic field. However, the higher MBSs may be sensitive to the polarization of 
the photon field. Moreover, comparing Figs. [1] and [2], we see that the energy dispersion 
is much more sensitive to the x-polarized photon field than that of ?/-polarization. This 
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Figure 1. (Color online) The many-body energy spectra for the interacting-electron 
and the quantized-photon modes versus the electron-photon coupling strength £ c for 
the case of TEon mode (the electric component with x-polarization). Shown are the 
pure photon states with no electron (green), the MBSs without the A 2 term (red), and 
the MBSs including the A 2 term (blue). The majority of the states contain a single 
electron. Two-electron states are present for energy > 2.0 meV. Other parameters are 
B = 0.1 T, W, a = 1.0 meV, and fiw = 0.4 meV. 

is because the single photon energy fkv is comparable to the characteristic energy of 
the lowest states in the nanostructure for electron motion in the ^-direction. In other 
words, it is caused by the anisotropy of the selected geometry. In addition, for the 
case of x-polarization the two states at around MBS energy Embs ~ 1.0 meV cross at 
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Figure 2. (Color online) The many-body energy spectra for the interacting-electron 
and the quantizcd-photon modes versus the electron-photon coupling strength £ c for 
the case of TEioi mode (the electric component with y-polarization). Shown are pure 
photon states with no electron (green), the MBSs without the A 2 term (red), and 
the MBSs including the A 2 term (blue). The majority of the states contain a single 
electron. Two-electron states are present for energy > 2.0 meV. Other parameters are 
B = 0.1 T, Ml a = 1.0 meV, and fiw = 0.4 meV. 

£ c pa 0.16 meV, but not in the case of y-polarization. 

We will now analyze the validity of the TLS approximation, as well as the even 
simpler JC model. The JC model is built from a TLS relying upon the assumptions 
of near resonance and weak coupling between the two systems that is described by the 
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Figure 3. (Color online) Comparison of the many-body energy spectra versus the 
coupling strength £ c for the case of TEon mode (a:-polarization). The energy spectra's 
are obtained by the TLS model with (blue) and without (green) the A 2 term and the 
JC model (red-dashed). The TLS model results are compared with the full numerical 
calculation results for the lowest active levels 1,2, and 3 (purple). Other parameters 
are B = 0.1 T, Ml = 1.0 meV, \g 12 \ = 0.701, and fkj = 0.185 meV. The inset shows 
the validity of the JC model in the weak coupling limit. 

following Hamiltonian in the second quantized form 

Hj C = -AEijO- z + huata + £jq (ct+ + a J) (a + a f ) , (22) 



Quantum magneto- electrodynamics of electrons embedded in a photon cavity 



11 



where AE^ = Ej — Ei denotes the energy difference between the electron states \i) and 
| j) which have been chosen as the relevant states for the TLS approximation. The ladder 
operators appropriate for a two-level approximation a± are defined by a± = ^(a x ±ia y ), 
where cr x ^ tZ are the Pauli matrices. Note that the energies of states |z) and \j) have 
been shifted to make them symmetric around the zero energy. 

The counter-rotating terms a + a' f and o_a in Eq. (1221) are usually omitted by taking 
the RWA to get an exactly solvable model. However, for our comparison we will keep 
the counter rotating terms and solve Eq. ( 1221) numerically using the Fock space basis 
{\k) <g) |M p h)}, where k G {i,j}- Comparison of the JC model with and without 
the counter rotating terms has been investigated [321 E31 EI], however it should be 
reexamined and compared with a system where realistic effects are included, such as 
those stemming from the non-trivial geometry of the nanostructure and an external 
magnetic-field. 

For a TLS with one electron, zero magnetic field and ignoring the S 2 term in Eq. 
(TT8l) . the total Hamiltonian in Eq. ([[]) reduces to a JC-like Hamiltonian with a coupling 
strength £ JC associated with £ c according to 

£jC = \9ti\£c, (23) 

in which the dimensionless coupling constant gij can be calculated using Eq. ( TTTj) . For 
a non-zero magnetic field Eq. (|23|) is not exact, however as will be shown later, a low 
magnetic field B — 0.1 T has minimal effects. 

Below we will assume that the detuning 5 is one percent of the energy spacing of 
the two active states AE^ giving the single-photon energy ftw = l.OlAEij. To label 
energy levels the notation Ejf is used. It refers to the energy of the state \k) ® |M p h) 
for 8 C = 0. 

Figure [3] shows the x-polarization many-body energy spectra as a function of the 
electron-photon coupling strength for the different models. We consider the lowest states 
|1) and |2) as the relevant states for the TLS. In the zero coupling limit £ c = 0, the 
ground state energy E\ = E® ~ 0.568 meV and the energy of the first excited state 
is E 2 = E\ — 0.751 meV so that the energy level spacing is AE U = 0.183 meV. The 
detuning is small compared to the typical energy difference of the electron states so 
the third state is associated with one-photon absorption from state |1) with energy 
E 3 = E\ = E® + ftw = 0.753 meV. As expected, the JC results almost coincide with our 
TLS results not including the A 2 term. The difference (not visible in Fig. |3]) between the 
two curves is due to effects of the external magnetic field. When the A 2 term is included, 
the energy spectrum manifests a blue-shift, and the energy- level correction SE may be 
larger than 0.02 meV in the strong coupling regime, or 5E/AE > 10%. A weaker red 
shift correction is observed when the higher MBSs are involved in the electron-photon 
coupling. 

In the weak coupling regime S c < O.lhui ~ 0.02 meV, the JC model is approximately 
valid. When the coupling strength is increased to £ c ~ hw ~ 0.2 meV, the ground state 
energy calculated by the TLS model is still valid. However, the energy of the excited 
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states becomes inaccurate, indicating that the simplified TLS model is no longer a 
good approximation in the strong coupling regime even though the diamagnetic vector 
potential A 2 is included. When the coupling strength £ c is increased, both the JC 
model and the TLS without the A 2 term predict a decreasing ground state, however by 
including the A 2 term within the TLS model the energy increases, in better agreement 
with our full numerical calculation. 

In Fig. HI we compare the many-body energy spectra as a function of the electron- 
photon coupling strength £ c when the electronic system is embedded in a TEjxu mode 
(y-polarization). Attributed to selection rules of the transverse parabolic confinement, 
we select the active states |1) and |5) to compare with the TLS approximation. In the 
zero coupling limit £ c = 0, we consider the ground state energy ~ 0.568 meV 
and the excited state energy El ~ 1.585 meV so that the energy level spacing 
A.E15 = 1.017 meV. In addition, we assume the detuning 5 15 = 0.01 x AE15 such 
that the single-photon energy fvjj = AE 15 + 5 15 = 1.027 meV. Moreover, we see 
that the state associated with a ground state electron absorbing one photon is around 
E\ = El + twj = 1.596 meV. 

Figure H] displays energy spectra calculated using the different models. As with the 
x-polarization, the energy spectrum obtained by the JC model almost coincides with the 
TLS result without the A 2 term. The difference (not visible in Fig. H]) is due to effects 
of the external magnetic field. When the A 2 term is included, the energy spectrum is 
blue-shifted in the strong coupling regime. When the higher MBSs are involved in the 
electron-photon coupling (full model) there is good agreement with the TLS including 
the A 2 term until inactive states (not included in the two-level approximation) start 
to have influence, such as the energy crossing at £ c ~ 0.13 meV and anti-crossing at 
£ c ~ 0.17 meV shown in Fig. HJ In the weak coupling regime £ c < O.lhuu ~ 0.1 meV, the 
JC model is approximately valid. The ultrastrong-coupling regime £ c > hui ~ 1.0 meV 
is not shown in this figure. 

The effects of the photon field on the charge distribution is illustrated in Fig. [5j 
where the charge distribution (Q(r)) of the third MBS is plotted for x-polarization for 
£ c = 0, 0.2, 0.3 and 0.4 meV. There is an energy crossing between the third and fourth 
state at £ c ~ 0.16 (see Fig. [1]) so it is important to note that for S c > 0.16, the third state 
refers to the fourth state counting from the bottom in Fig. [TJ In Fig. |5^a), the energy of 
the MBS is Eq ~ 0.96 meV. There is no coupling between the photons and electrons so 
the charge density is identical to that of the ground state labeled by the energy E®. For 
£ c = 0.2 meV= fiui/2, the charge distribution is stretched in the x-direction. This trend 
continues with increasing coupling strength. For £ c = 0.3 meV the charge distribution 
starts to separate into two peaks and at £ c = 0.4 meV the two peaks are completely 
separated. In other words, in the ultrastrong coupling regime £ c = froj, a clear dipole- 
like charge distribution profile is observed. For the ^/-polarization, not shown here, the 
polarization of the charge is much smaller due to the large value of the confinement 
energy HQq = 1-0 meV compared to the photon energy fiw = 0.4 meV. The system is 
anisotropic at the energy scale we employ here. 
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Figure 4. (Color online) Comparison of the many-body energy spectra versus the 
coupling strength £ c for the case of TEioi mode (y-polarization). These energy states 
are obtained by TLS model including the A 2 term (blue), not including the A 2 term 
(green) and the JC approximation without magnetic field (red-dashed). The TLS 
model results are compared with the full numerical calculation for the compared lowest 
active levels 1, 5, and 6 (purple) as well as inactive levels (gray). Other parameters 
are B = 0.1 T, M7 = 1.0 meV, |pi, 5 | = 0.290, and hu = 1.027 meV. The inset shows 
the validity of the JC model in the weak coupling limit. 

To summarize this section, we remind the reader that the JC model exhibits 
energy spectrum with some levels decreasing when the electron-photon coupling is 
increased [321 [33]. The energies may be negative if the electron-photon coupling is 
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Figure 5. (Color online) Charge distribution (in units of — e) of the third many- 
body state (see text for definition) in the case of s-polarization with clcctron-photon 
coupling strength £ c = (a) 0.0, (b) 0.2, (c) 0.3, and (d) 0.4 meV. The other parameters 
are the same with Fig. □ B = 0.1 T, Ml = 1.0 meV, and hw = 0.4 meV. 

very strong [32J . Comparing with the full numerical calculation, it is thus unambiguous 
that the diamagnetic A 2 contribution, as well as higher energy electron states have to be 
included in the ultrastrong-coupling regime. We have done a simple literature search for 
"ultrastrong" and "circuit-QED" . One has to keep in mind that authors have not always 
used the same definition for strong or ultrastrong, but none found in our search results 
included the A 2 term in their model nor went beyond the two-level approximation. 

4. Concluding Remarks 

We have performed a numerical calculation of a microscopic model describing a hybrid 
structure consisting of an electronic nanostructure embedded in a cavity resonator. 
We have demonstrated strong coupling features of Coulomb interacting electrons and 
photons in a nanostructure embedded in a cavity resonator in an external magnetic 
field. The two-dimensional electronic nanostructure is parabolically confined in the y- 
direction and hard-wall confined in the x-direction that is embedded in a rectangular 
photon cavity with a TE-mode electromagnetic field that may be either x- or y-polarized. 
We have found that the many-body energy spectrum is more sensitive to the photon 
field with x-polarized electric component than that with y-polarization for the selected 
geometry. The system is anisotropic in the energy range explored. 

We have established that the diamagnetic A 2 term in the Hamiltonian may provide 
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a blue-shift correction to the energy spectrum. However, when higher many-body 
states are included beyond a two-level approximation, the results of the full numerical 
calculation exhibit a smaller red-shift correction. This implies that the two lowest levels 
become more stable when the higher energy levels are included in the electron-photon 
coupled system. When the A 2 term is not included the energy spectrum is lowered when 
the coupling strength is increased, but the opposite trend is found when the A 2 term is 
included in the calculation. 

The widely employed two-level system approximation has been reexamined 
comparing to results of our full numerical calculation model. Qualitative difference of 
the energy spectrum between the JC model and the full numerical calculation is found 
in the strong-coupling regime. The JC-model includes no information about the charge 
distribution of the system. A strong cavity photon field can cause large polarization of 
the charge distribution, an effect seen in Fig. [5j The reason for the high number of single- 
electron states needed in our full calculation is exactly this large polarizing effect of the 
photon field. QED modeling of a circuit element on the nanoscale in the ultrastrong 
coupling regime requires approximations beyond the JC-model or more general two level 
models. 

In summary, we have presented a model adequate for accurate numerical calculation 
for the electron-photon coupled energy spectrum that is essential and will be utilized to 
explore time- dependent transport of electrons through a photon cavity in a forthcoming 
publication [31 J . 
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